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Abstract. Once distributed all over the Philippines, the endemic Philippine Crocodile (Crocodylus mindorensis) is nowadays 
threatened with extinction. It is estimated that less than 140 mature individuals live in the wild. Human activities like fish- 
ing and poaching, as well as land-use change and habitat conversion cause a continuing threat to the remaining populations. 
Therefore, designated protected areas (PAs) were evaluated with species distribution models (SDMs) and also to see if most 
suitable areas are covered by PAs in order to improve future conservation efforts. For this purpose, the existing IUCN-re- 
serves were analysed for potential habitat suitability (combining bioclimatic and remote sensing variables), wetland occur- 
rences and the human footprint index by using MaxEnt and QGIS. Based on species records, our final SDM showed high 
performance and revealed the climatically most suitable areas for the species, which were mostly on Luzon and Mindanao. 
However, only small parts of the climatically suitable wetlands are currently covered by reserves (0.3-46.3%). In addition, 
none of the species’ records was located within a PA. The anthropogenic pressures in the reserves measured by human foot- 
print index (considering eight variables i.e. ‘population density, ‘navigable waterways, ‘crop lands’ and ‘roads’) were diverse 
and varied between a low and moderate level. Most of the records were found in areas with a moderate human footprint. 
Considering the three criteria, ‘Lake Lanao Watershed Reservation, ‘Angat Watershed Forest Reserve District (Metro Wa- 
ter District), ‘Northern Sierra Madre Natural Park, “Talaytay Protected Landscape’ and ‘Agusan Marsh Wildlife Sanctuary’ 
revealed to be the most suitable conservation areas for C. mindorensis, whereas suitable areas outside PAs are highly recom- 
mended for further surveys. We recommend to declare Ligawasan Marsh, Mindanao as a PA as this area harbours a large 
population of C. mindorensis. The declaration of more climatically suitable areas with low level of human footprint to PAs is 
a necessary step for the long-term conservation of this endemic crocodile species. The current network of existing PAs needs 
improvement in order to provide well-suited and long-term protection for C. mindorensis. More surveys are also necessary 
to find hidden, so far overlooked populations and to assess C. mindorensis tolerance level for human impacts. 
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Introduction 


The Philippines are one of the 20 global megadiverse coun- 
tries and a major biodiversity hotspot in Southeast Asia 
due to its isolated location and diverse topography (CBD 
2018, PERIA 2014, UNDP 2021, VON RINTELEN et al. 2017). 
These megadiverse countries are home to about 70-80% of 


the plant and animal species on the planet, of which more 
than 20,000 are endemic (AMBAL et al. 2012, FPE 2013). 
Almost half of the terrestrial animals occurring in the Phil- 
ippines are also endemic to the country and in the case of 
reptiles, ~ 70% (244 of the 352 known species in 2017) of 
native species are endemic (PSA 2019). According to the 
categories of the International Union for Conservation of 
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Nature’s (IUCN) Red List of Threatened Species, 652 native 
species of animals and plants are classified as “Vulnerable; 
542 as ‘Endangered, 309 as ‘Critically Endangered’ and 15 of 
those (all members of the fish genus Barbodes) are already 
extinct (IUCN 20214). 

The Philippine Crocodile Crocodylus mindorensis 
Schmidt, 1935 has been classified as ‘Critically Endangered’ 
on the IUCN Red List since 1996, while the latest assess- 
ment was in 2012 (VAN WEERD 2016). Unlike the Indo-Pa- 
cific Crocodile, Crocodylus porosus Schneider, 1801, also 
known as Saltwater Crocodile, which is native to the Phil- 
ippines but extends as far as South Asia, Southeast Asia 
and Australia (UETZ 2021), the Philippine Crocodile is en- 
demic to the Philippines (VAN DE VEN et al. 2017). Origi- 
nally, C. mindorensis was distributed over almost the entire 
archipelago and inhabited islands of Mindoro, Masbate, 
Samar, Negros, Busuanga, Luzon and Mindanao (Ross 
1982, Ross & ALCALA 1983). Since the late 1990s and in the 
2000s, more surveys have been conducted, but informa- 
tion on the actual distributional range of the crocodile re- 
mains scarce. Although it is protected by law since 2001, 
there are only about 92-137 estimated mature individuals 
left in the wild to date (VAN WEERD 2016), and populations 
are estimated to decline. These are highly fragmented pop- 
ulations in ‘Dalupiri Island, ‘Northern Luzor’ and ‘Liga- 
wasan Marsh on Mindanao (MANALO et al. 2013, MANALO 
et al. 2015, VAN WEERD 2016). 

The habitats of the relatively small Philippine Crocodile 
are wetlands with freshwater occurrences such as creeks, 
ponds, man-made water reservoirs, mangrove areas and 
marshes, but also fast-flowing rivers with caves made from 
limestone cliffs (VAN DE VEN et al. 2017). These caves are 
used as hiding places just as the ones in sandy and clay 
river banks. Similar behaviour has also been documented 
with the introduced population of the species in Paghun- 
gawan Marsh, Siargao Island which was part of the govern- 
ment’s effort to repopulate the species in the wild (BINADAY 
et al. 2020). The species’ altitudinal range extends from the 
favoured inland wetlands up to 850 m (sea level vs. Cordil- 
lera Mountains on Luzon) (MANALO 2007). 

The Philippines face several environmental problems 
like deforestation and forest degradation, water pollution, 
poaching and illicit wildlife trade (USAID 2021). The na- 
tional desire for more sustainability, environmental pro- 
tection and species conservation is often contrasted by the 
poor income situation of local communities (ADAMs et al. 
2004, JAISANKAR et al. 2018). Furthermore, large parts of 
protected areas (PAs) overlap with the ancestral domains 
(PEREZ 2018). Local communities living close to or even 
within these areas rely on the local resources and will be 
socio-economically harmed by strict environmental reg- 
ulations, especially if there are no alternative livelihoods 
(ApaAMs et al. 2004). In fact, some of the greatest threats 
for the Critically Endangered C. mindorensis is the use of 
its natural habitat by rural people, as well as habitat de- 
struction. In addition, the crocodiles are often persecuted 
and their nests destroyed or plundered by humans. Fishing 
is also considered a danger to these crocodiles as they are 
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likely caught in fishing nets as bycatch (AKMAD & POMARES 
2008, VAN WEERD 2016). 

As early as 1992, the Philippine government committed 
itself to the international goals of the Convention on Bio- 
logical Diversity (CBD). This resulted in several national 
environmental laws such as the National Integrated Pro- 
tected Area System (NIPAS) Act of 1992 and the Wildlife 
Resources Conservation and Protection Act (2001), which 
are intended to protect the country’s natural resources in 
the long term (DENR-BMB 2021). Currently, 248 areas 
have been recognised by the Philippine government as PAs, 
covering a cumulative area of 7.8 million ha (DENR-BMB 
2020). However, it is currently unknown if these PAs pro- 
vide climatically suitable habitats for the Philippine Croco- 
dile. Species distribution models (SDMs) have been widely 
used and proven to be very useful in habitat analyses of 
other species and are used for prioritisation in conserva- 
tion planning (Brnapay et al. 2020, Fots et al. 2018, IHLOW 
et al. 2015, RODDER et al. 2010, TAN et al. 2022, TSUYAMA 
et al. 2015). SDMs attempt to predict potentially environ- 
mentally suitable habitat by linking documented presence 
records of species to environmental variables and spatial 
characteristics such as human footprint and availability of 
surface water based on the species’ ecological niche. In this 
study, it was investigated whether the existing PAs (1) pro- 
vide suitable wetland habitats and (2) are climatically suit- 
able for the Philippine Crocodile. Furthermore, we (3) in- 
cluded anthropogenic impact measured as human foot- 
print index to identify these PAs where low anthropogenic 
pressure occurs. 


Methods 


For the evaluation of suitable PAs for the species, the 248 
current PAs, availability of wetland areas, the climatic suit- 
ability and anthropogenic pressure were considered. Since 
the first two criteria are decisive for the basic survival of 
the ectothermic species, the final ranking involved three 
steps. In the first step, the wetlands were evaluated. In a 
second step, the climatically suitable areas of the remain- 
ing 117 sites were identified. Anthropogenic pressure to 
the remaining 114 sites was assessed in the third step. A 
final ranking was calculated based on the combined pro- 
portions of suitable wetland area and climate suitability 
(‘wet-sdm-ranking’) and anthropogenic pressures (‘hfp- 
ranking’), which were multiplied with each other. Below 
we describe our workflow in detail. 


PAs and wetlands data 


As it is easier to implement in situ conservation measures 
of C. mindorensis in areas that are subject to minor an- 
thropogenic influence, the coverage of the species range 
with PAs was assessed in addition to the assessment of 
available potential habitats. PAs are defined as geographi- 
cal areas and classified by IUCN standards (categories I- 
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VI) to achieve the long-term conservation of nature and 
the corresponding ecosystems (IUCN 2021b). The World 
Database of Protected Areas (WDPA) polygon shapefiles 
were obtained from UN Environment Programme World 
Conservation Monitoring Centre (UNEP-WCMC; UNEP- 
WCMC & IUCN 2021). 

The importance of wetlands was evaluated based on as- 
sessed tropical and subtropical wetlands with a resolution 
of 232 meters downloaded from the Global Wetlands Map 
(https://www2.cifor.org/global-wetlands/) (MANALO et al. 
2018, VAN WEERD & VAN DER PLOEG 2012). GUMBRICHT et 
al. (2017) developed a mapping method combining differ- 
ent data sources and methods, and classified wetlands into 
three key biophysical attributes: ‘long-term water supply, 
‘annually or seasonally water-logged soils’ and ‘a geomor- 
phological position where water can be supplied and re- 
tained’. Seven categories were selected for C. mindorensis: 
‘open water, ‘mangrove, ‘riverine, ‘floodplains, ‘marshes, 
‘swamps and ‘fens. Furthermore, an additional category 
was added by importing a high-resolution water layer from 
the Global Runoff Data Centre (GRDC) to assess the suit- 
ability of river networks for the crocodile (GRDC 2020). 


Species records and predictor variables 


Species occurrence records were collected by JWB and RM 
between 2003 and 2021. The dataset was examined for out- 
liers in QGIS, ver. 3.16.3 with GRASS 7.8.5 (QGIS.org 2021), 
but not corrected for potential spatial autocorrelation due 
to the few occurrence records. Habitat suitability was pre- 
dicted using SDMs based on a combination of 46 environ- 
mental variables ( Supplementary Table S2; Corp & R6p- 
DER 2011). The 19 bioclimatic variables were obtained from 
Worldclim database, ver. 1.4 and contain interpolated el- 
ements from different climate conditions collected over a 
period of 30 years (1960-1990) with a resolution of 30 arc 
seconds (HIJMANS et al. 2005). The remaining 27 environ- 
mental predictors were derived from Moderate Resolu- 
tion Imaging Spectroradiometer (MODIS) sensors of two 
NASA satellites. The spatial resolution of the pre-processed 
remote sensing variables amounts to 30 arc seconds and 
the temporal resolutions are 8-day averages (MOD11A2) 
and 16-day averages (MCD43B4) (Mu et al. 2007, SCHAR- 
LEMANN et al. 2008). Since SDMs are sensitive to multi- 
collinearity of predictors (DE MARCO & NÓBREGA 2018, 
Merow et al. 2013), we calculated Variance Inflation Factor 
(VIF) in R (ver. 4.0.3; ‘usdm package’) (Narmı et al. 2014, 
R Core Team 2020) to exclude highly correlated variables, 
when one of them exceeded the value of 10. The final vari- 
ables were temperature ranges, precipitation and isother- 
mality ( Supplementary Table S2). 


Species distribution modelling 


For SDM, MaxEnt was chosen as this machine-learning 
programme (PHILLIPS et al. 2006, PHILLIPS et al. 2016) is 


shown to be more reliable than other modelling tools es- 
pecially when dealing with small sample sizes (EL1TH et al. 
2006). As the historic distributional range of C. mindoren- 
sis covers large parts of the Philippines (UETZ 2021), the 
whole country was chosen as background area. In addition, 
the records used for SDM construction were reduced to one 
per grid cell to reduce sampling bias (PHILLIPs et al. 2009). 

Model fitting and selection followed the procedure de- 
scribed in GINAL et al. (2022) and is based on testing mul- 
tiple regularisation multipliers (0.5 to 2.5 in steps of 0.1, as 
well as 5 and 10) and feature classes (LP, LQ, LH, LT, LQP, 
LQH, LQT, LPH, LPT, LHT, LQPT, LQHT, LPHT, LQPHT; 
L = Linear, P = Product, Q = Quadratic, H = Hinge, T = 
Threshold)). MaxEnt’s raw output format was used for 
further processing and model selection, and the averaged 
AICc [corrected Akaike Information Criterion (WARREN & 
SEIFERT 2011)] and AUC (ELITH & GRAHAM 2009, LOBO et 
al. 2008, PHILLIPS & DUDÍK 2008) were calculated across 
ten replicates. Further, AUC was used as an evaluation of 
the model performance (ELITE et al. 2010). For AUC calcu- 
lation, the presence data were randomly divided for model 
training (80%) and testing (20%) using the bootstrap ap- 
proach. For model selection, the lowest average AICc and 
an AUC. „above 0.7 were used (PHILLIPS & DUDÍK 2008, 
WARREN & SEIFERT 2011). The final model was replicated 
100 times, again with an 80:20 split for training and test- 
ing. Finally, the average over the 100 replicates was calcu- 
lated and evaluated using a combination of AUC (ELITH 
et al. 2006) and True Skills Statistics (TSS) (ALLOUCHE et 
al. 2006, SHABANI et al. 2018). For the final model, cloglog 
format was used as output. Considering the limited num- 
ber of available occurrence records and the historical dis- 
tribution of the species, the ‘minimum training presence’- 
threshold was chosen for presence/absence. 

The above mentioned wetland shapefile was overlaid 
with the reclassified MaxEnt-output (settings ‘o - thresh- 
old = NA; = threshold = r) and then analysed together 
with PAs using the ‘zonal.histogram’-raster function in 
QGIS (QGIS.org 2021). The obtained numbers of grid cells 
per category were summed up per reserve in proportion to 
the total area of the reserve. 

To obtain sums and counts from the final model, the 
MaxEnt output was subjected to a second but separate 
classification in a first step (settings ‘o - threshold = NA). 
Then this reclassified MaxEnt output was rescaled within 
the range o-1 before it was analysed with the shapefile gen- 
erated in the ‘zonal.histogram-analysis’ using the ‘zonal. 
statistics’-raster function in QGIS. The generated data pro- 
vided information on how well the habitat is suited for 
the Philippine Crocodile in terms of climatic conditions 
(‘sum’). Furthermore, it was possible to calculate the area 
of suitable habitat within a PA using ‘count, which calcu- 
lates the number of grid cells of the suitable area. The sum 
values were ranked in descending order. Since the reso- 
lutions of the MaxEnt-map (~ 1000 m) and the wetlands 
map (232 m) differ, both rankings were multiplied and a 
new combined ranking was assigned (‘wet-sdm-ranking, 
ascending order). 
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Table 1. Results of the ten best MaxEnt models used for model 
selection, ranked by the mean AICc values and with informa- 
tion on the regularisation multipliers, feature classes, number of 
parameters, AICc, AUC,, ,, and AUC hse The final model used for 
the following processes is shown in bold. 


Regulari- Features nParameters AICc AUC AUC 


5 Train Test 

sation 

0.6 LPT 7.5 253.40 0.92 0.86 
0.9 L 6.5 256.08 0.90 0.83 
1.0 LP 6.5 260.48 0.93 0.89 
0.8 LP 6 265.08 0.88 0.89 
1.2 LPT 5.5 266.02 0.88 0.86 
1.1 LPT 6 266.06 0.86 0.84 
1.0 LT 7 266.30 0.90 0.87 
1.1 L 5.5 266.67 0.87 0.83 
1.0 L 6.5 266.94 0.88 0.83 
1.3 LT 5 267.03 0.90 0.85 


Anthropogenic pressure 


To assess the potential effect of anthropogenic pressure on 
the crocodiles, the 2018 release of human footprint was ob- 
tained from SEDAC (NASA Socioeconomic Data and Ap- 
plications Center) and added to our analyses. These maps 
comprise eight variables (i. e. ‘population density; ‘naviga- 
ble waterways; ‘crop lands’ and ‘roads’) to measure the di- 
rect and indirect human pressure (VENTER et al. 2018). The 
human footprint-ranking (hfp-ranking’) was computed 
based on means of the “zonal.statistics’-raster function of 
the combined map of “wildareas v3 2009 human footprint” 
and previously mentioned PAs shapefile. 


Results 
Species distribution modelling 


For model fitting, MaxEnt computed 3450 models in to- 
tal (23 regularisation multipliers x 15 feature class com- 
binations x 10 replicates) of which the ten best perform- 
ing models were ranked according to the lowest average 
AICc (Table 1). The ten models revealed high AUC values 
[AU C. ain 0.86 - 0.93, AUC. _, 0.83 — 0.89]. The values of the 
final model, which was replicated 100 times, were: regu- 
larisation multiplier 0.6, feature classes LPT, AUC,,,, 0.92, 
AUC,,,, 0.86 and TSS 0.45 + 0.18. ‘Mean diurnal range of 
temperature’ had the highest contribution to the final SDM 
(36.4%), followed by ‘isothermality’ (16.8%), ‘seasonality’ 
(8.5%), ‘precipitation of coldest quarter’ (8.4%) and ‘annual 
range of NDVI (7.9%). The remaining variables contribut- 
ed only less to the model performance (Table 2). 


Availability of wetlands, climatic suitability and 
anthropogenic pressure 


Only 57 of the 248 national PAs are currently designated by 
the IUCN. The areas of the 248 PAs strongly differed and 
ranged from 0.04 km? (‘HinuluganTaktak Protected Land- 
scape, Luzon) to 10,881.81 km? (‘Palawan Game Refuge and 
Bird Sanctuary; Luzon). In the first step of ranking, 131 PAs 
were excluded from further evaluation due to the lack of 
habitat availability (Supplementary Material Si). In the 
second step another three areas were excluded as they did 
not provide suitable climatic conditions. The remaining 114 
PAs were included in the final ranking. 

Considering the results of the ‘sdm-ranking, the PAs 
‘Northern Sierra Madre Natural Park’ (North Luzon), 


Table 2. MaxEnt variable contribution of the final species distribution model for Crocodylus mindorensis. 


Variable Abbreviation Derived variable Variable contribution [%] 
V39 ED15078_bio2 Mean Diurnal Range of Temperature 36.4 
V40 ED15078_bio3 Isothermality (Bio2/Bio7) (*100) 16.8 
V41 ED15078_bio4 Seasonality 8.5 
V19 bio_19 Precipitation of Coldest Quarter 8.4 
V26 ED1514_bio7 Annual Range of NDVI 7.9 
V37 ED1515_bioll Mean EVI of Coldest Quarter 5.2 
V14 bio_14 Precipitation of Driest Month 4.2 
V18 bio_18 Precipitation of Warmest Quarter 3.9 
V35 ED1515_bio7 Annual Range of EVI 3.8 
V27 ED1514_biol0 Mean NDVI of Warmest Quarter 2.9 
V13 bio_13 Precipitation of Wettest Month 1.9 
V43 ED15078_bio6 Min Temperature of Coldest Month 0.1 
v31 ED1515_bio3 Isothermaility (BIO2/BIO7) (*100) of EVI 0.0 
V30 ED1515_bio2 Mean Diurnal Range of EVI 0.0 
V22 ED1514_bio3 Isothermaility (BIO2/BIO7) (*100) of NDVI 0.0 
V21 ED1514_bio2 Mean Diurnal Range of NDVI 0.0 
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Table 3. Top 5 reserves suitable for the Philippine Crocodile: reserve name, reserve category (assignment according to IUCN), reserve 
area, climatically suitable area [relative to reserve area in %], climatically suitable wetland area [relative to reserve area in %], and 
ranks according to the wet-sdm-ranking, hfp-ranking, and final-ranking. 


Name Reserve category IUCN reserve  climatically wetland wet- hfp- final- 
area suitable area area[km?] sdm- ranking ranking 
[km?] [km?] ranking 

Lake Lanao Watershed Watershed not 1712.93 946.87 113.82 2 34 1 

Reservation Reservation assigned [55.3%] [6.6%] 

Angat Watershed Forest Reserve Watershed Forest not 545.74 191.96 2.83 38 2 2 

District (Metro Water District) Reserve assigned [35.2%] [0.5%] 

Northern Sierra Madre Natural Natural Park II 3569.69 1664.39 96.22 1 86 3 

Park [46.6%] [2.7%] 

Talaytay Protected Landscape Protected V 35.98 1.44 0.09 100 1 4 

Landscape [4.0%] [0.3%] 
Agusan Marsh Wildlife Wildlife IV 409.41 247.32 189.43 3 37 5 
Sanctuary Sanctuary [60.4%] [46.3%] 


‘Lake Lanao Watershed Reservation (West Mindanao), 
‘Quirino Protected Landscape (Luzon), Allah Valley Wa- 
tershed Forest Reserve (South Mindanao) and ‘Upper 
Agno River Basin Resource Reserve (Luzon) revealed the 
highest scores with climatically suitable areas between 
549.29 and 1,664.39 km’ (Table 3, Supplementary Material 
S1). 

According to the ‘wet-sdm ranking, the top five re- 
serves with the highest scores were ‘Northern Sierra 
Madre Natural Park (North Luzon), ‘Lake Lanao Water- 
shed Reservation’ (West Mindanao), Agusan Marsh Wild- 
life Sanctuary’ (East Mindanao), ‘Mindoro Island’s Man- 
grove Swamp Forest Reserves as per Presidential Procla- 
mation 2152’ (South Luzon) and ‘Allah Valley Watershed 
Forest Reserve’ (South Mindanao). This coincided with 
the distribution of the species records, which were also 
identified on the Northern portion of Luzon Island and 
Mindanao Island (Figs 1 and 2). However, none of the spe- 
cies’ occurrence records laid inside the top five reserves. 
Only two records were located on the edge or close to a 
reserve (‘Northern Sierra Madre Natural Park). Follow- 
ing the IUCN categories, three PAs were not assigned to 
any IUCN category (‘Lake Lanao Watershed Reservation, 
‘Mindoro Island’s Mangrove Swamp Forest Reserves’ and 
‘Allah Valley Watershed Forest Reserve’), while Agusan 
Marsh Wildlife Sanctuary’ belongs to category IV ‘habi- 
tat/species management area, and ‘Northern Sierra Madre 
Natural Park’ is classified as ‘national park (category II). 
‘Lake Lanao Watershed Reservation and ‘Allah Valley Wa- 
tershed Forest Reserve are two reserves proclaimed by the 
national government through Presidential Proclamations 
No. 871 and 2455, respectively. Governance and manage- 
ment of these reserves are also covered by the NIPAS Act 
of 1992. Meanwhile, the Presidential Proclamation 2152 de- 
clares several mangrove areas throughout the country as 
‘Mangrove Swamp Forest Reserves; this includes the man- 


groves areas of Mindoro Island mentioned in this study. 
Our analysis showed that there are generally few areas 
with low anthropogenic pressure except for the mountain 
ranges on Luzon (Fig. 3). Not surprisingly, high human ac- 
tivity was found around the capital Manila. On the main 
island of Palawan, the human footprint was low, whereas 
in the Visayas, except for ‘Samar Island’ there were only a 
few contiguous areas with low human footprint. Minda- 
nao, meanwhile, has a very homogeneous pattern distrib- 
uting between high and low anthropogenic pressure. The 
reserves with the lowest human footprint were “Talaytay 
Protected Landscape’ (Central Luzon, IUCN category V 
= ‘protected landscape/seascape’), ‘Angat Watershed For- 
est Reserve District (Metro Water District)’ (Luzon, not 
assigned), ‘Mt. Mantalingahan Protected Landscape’ (Lu- 
zon, IUCN category V), Amro River Protected Landscape’ 
(Central Luzon, IUCN category V) and “Mt. Pulag Pro- 
tected Landscape’ (Luzon, not assigned) (Supplementary 
Material S1). 

Considering the availability of wetlands, the climat- 
ic suitability and the anthropogenic pressure, the final 
ranking revealed ‘Lake Lanao Watershed Reservation 
(West Mindanao, not assigned), ‘Angat Watershed For- 
est Reserve District (Metro Water District) (Luzon, not 
assigned), “Northern Sierra Madre Natural Park’ (North 
Luzon, IUCN category II), “Talaytay Protected Landscape’ 
(Central Luzon, IUCN category V) and ‘Agusan Marsh 
Wildlife Sanctuary’ (East Mindanao, IUCN category IV) 
as most suitable reserves for the Philippine Crocodile. The 
reserves covered a total area of between 35.98 km” and 
3,569.69 km?, with climatically suitable areas of between 
1.44 km? and 1,664.39 km’, and finally consisted of be- 
tween 0.09 km? and 189.43 km? climatically suitable wet- 
land habitats. The anthropogenic pressure strongly varied 
between low (ranks 1, 2, 34 and 37 in the ‘hfp-ranking’) 
and moderately high (ranks 86; Tables 3 and S1). 
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Discussion 
General results 


The assessed PAs strongly differ in size, IUCN/reserve cat- 
egory, climatic suitability, habitat availability and anthro- 
pogenic pressure, and therefore a trade-off is necessary 
to identify the most suitable PAs to protect the Philippine 
Crocodile in situ. Our SDMs indicated an average daily 
temperature range (36.4%), temperature fluctuations dur- 
ing the course of the day and year (16.8%), and seasonality 
(8.5%) to be most relevant climatic parameters for the Phil- 
ippine Crocodile (Table 2). Ambient temperature influenc- 
es the nest temperatures for reproduction and long-term 
maintenance for the species (AKMAD & POMARES 2008). 


The amount of precipitation in the coldest quarter was 
also found important for the crocodiles (Table 2). Specific 
weather or climatic events such as floods can cause mor- 
tality to juveniles when they are still vulnerable (VAN DE 
VEN et al. 2009, VAN DE VEN et al. 2017). Furthermore, pro- 
longed dry periods limit food resources and are therefore a 
threat for all age groups (MazzoTTI! et al. 2009). 


Assessment of top five PAs and recommendations to 
establish new reserves 


‘Lake Lanao Watershed Reservation (not assigned) is a wa- 
tershed reservation located in the province of Lanao del 


@ Records 
O Protected areas 
Wetland 


Figure 1. Map of the Philippines including species records of the Philippine Crocodile (Crocodylus mindorensis), wetland habitats, 


and national PAs. 
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Sur of the Bangsamoro Autonomous Region in Muslim 
Mindanao (BARMM) on the island of Mindanao. With a 
total reserve area of 1,712.93 km’, covering 113.82 km? of 
climatically suitable wetland habitats, it is the second larg- 
est of the top five. Lake Lanao is Mindanao’ largest lake 
(36.300 ha) and has five watersheds with rivers and major 
tributaries stretching over a total length of 431 km (DENR 
2023). Its wetlands border the lake to the east for the most 
part and mainly consist of general marshes (39%), swamps 
(26%), fens (16%) and rivers (14%). The hydropower plant 
built along the lake and the Agus River is responsible for 
a significant contribution to Mindanaos electricity sup- 
ply (70%) (DENR 2023), which also reflected by high an- 
thropogenic pressure. The PA is also recognised as key bio- 
diversity area (KBA) by IUCN and is therefore of crucial 
global importance. The lake is home to 18 endemic fresh- 
water fish and supports a large number of waterfowl and 
other birds such as Halcyon chloris (White Collared King- 
fisher) (DENR 2023). Moreover, a healthy population of 


C. mindorensis is inferred to be thriving in the rivers of 
Miundas, Maladi and Matling in Lanao del Sur with a re- 
cent discovery of individuals in 2019 and affirmation of its 
presence by the local community (MANALO et al. 2019). 
The headwaters of these three rivers are located in the vi- 
cinity of Lake Lanao. 

The ‘Angat Watershed Forest Reserve District (Metro 
Water District)’ (not assigned to IUCN categories) pro- 
tects the watershed of the Southern Sierra Madre north 
of Manila, where surface water flows into the Angat River 
and its tributaries. The rivers hold a proportion of 74% of 
the total wetlands in the PA. The reserve covers an area of 
545.74 km? in the eastern part of Bulacan Province and the 
northern portion of Rizal Province at elevations between 
490 to 1,206 metres a.s.l. The PA extends to the provinces 
of Nueva Ecija and Quezon and is centred on an artifi- 
cial lake created by the Angat Dam which, together with 
the Ipo Dam (7.5 km downstream), supplies the majority 
of Metro Manila’s water requirements. Despite this fact, 
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Figure 2. Suitable climatic space for the Philippine Crocodile (Crocodylus mindorensis) according to the results of our species distri- 
bution model as well as species records and national PAs. The five best suitable PAs according to our final ranking are highlighted. 
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the PA has a homogenous low human footprint. The wa- 
tershed is a popular birdwatching site and a biodiversity 
hotspot, containing most of the remaining closed cano- 
py forests in Central Luzon. A herpetofaunal survey con- 
ducted by Mctzop et al. (2011) documented 19 frog, 22 
lizard, two turtle, and 20 snake species within the PA, but 
survey efforts were mostly focused on low elevation sites 
(200-600 m a.s.l.). However, the presence of C. mindoren- 
sis in this PA is still unknown and needs further surveys. 
Although the area is highly suitable for the crocodile, large 
dams may restrict movement of the animals and the sepa- 
ration of populations would have a detrimental effect on 
the long-term conservation of the species (MCALLISTER 
et al. 2001) or would require management or assisted mi- 
gration. Nevertheless, the Philippine Crocodile has been 
documented to be able to climb steep slopes (BINADAy et 
al. 2020) and studies are yet to be done on whether such 
infrastructures will have a significant impact on the spe- 
cies’ population. 


‘Northern Sierra Madre Natural Park, a large national 
park, covers 3,569.69 km? of the north-east coast of Lu- 
zon. About 21.02% of the park is highly suitable. The wet- 
land area, covering 17 km”, is mainly consisting of rivers, 
swamps and marshes. The areas of the Sierra Madre Moun- 
tains on the east coast and the Cordillera Mountains on 
the middle-west side are exposed to little or no anthropo- 
genic pressure. Only the settlements and human activities 
along the branches of the Palanan River possess medium 
to high human footprint. In the east of San Mariano, there 
are already several crocodile sanctuaries for this species 
(MANALO et al. 2018, VAN DE VEN et al. 2017). However, 
the sanctuaries are currently located outside the PA where 
three of the species occurrences were found (Fig. 1). We 
highly recommend the extension of the designated reserve 
to cover areas surrounding the sanctuaries and especially 
the area north along the foothills of the Dicatian River as 
a potential reserve, where there is high climatic suitability 
and low human footprint (Figs 2 and 3). This park contains 
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Figure 3. Anthropogeneic pressure in the Philippines measured as human footprint index as well as species records of the Philippine 
Crocodile (Crocodylus mindorensis) and national PAs. The five best suitable PAs according to our final ranking are highlighted. 
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the highest species richness of the Philippines, a wide va- 
riety of habitats (DENR 2015) and is home to many indig- 
enous people (CEPF 2001). Among the Philippine endemic 
species, about 30% of all bird species and 62% of all mam- 
mal species are found here. It is also home to 35 threatened 
species (VAN DER PLOEG et al. 2011), including C. mindo- 
rensis, Philippine Eagle (Pithecophaga jefferyi), Isabela Ori- 
ole (Oriolus isabellae) and Sierra Madre Forest Monitor 
(Varanus bitatawa) (DENR 2015). A herpetofaunal survey 
was conducted by Brown et al. (2013) in the PA and docu- 
mented a total of 101 species of amphibians and reptiles, 
including the two species of crocodiles in the Philippines - 
Crocodylus mindorensis and C. porosus. Although much of 
the PA is covered by forest, it is reducing by about 1,400 ha 
per year (DENR 2015). VAN DER PLOEG et al. (2011) esti- 
mated that between 20,000-35,000 m° disappear from the 
national park each year due to illegal timber logging, but 
little action has been taken against this so far. 

“‘Talaytay Protected Landscape’ is located in northern 
Aurora (province) and covers the Talaytay River watershed 
in the Sierra Madre range of the island of Luzon. The PA 
comprises an area of 35.98 km’, making it the smallest of 
the top five PAs, but like ‘Lake Lanao Watershed Reserva- 
tion’ it is a KBA. This Protected Landscape (IUCN cate- 
gory V) stretches from the rugged interior including the 
source region of the Talaytay River to its mouth at the low- 
lands of the municipality of Dinalungan. Some important 
bird and mammal species are native to this PA, such as 
Penelopides panini (Tarictic Hornbill) or Macaca fascicula- 
ris philippensis (Philippine Long-Tailed Macaque). The 
wetlands there consist exclusively of rivers, and anthropo- 
genic pressure is low. ‘Northern Sierra Madre Natural Park’ 
and “Talaytay Protected Landscape are PAs that are close to 
the coast or include parts of it. These habitats might also be 
suitable for the much larger C. porosus, which might out- 
compete the Philippine Crocodile for food, nesting sites or 
basking sites. Furthermore, hybrids between both species 
are known from captivity, which should be considered for 
conservation actions in these areas. 

Almost 60% of Agusan area are climatally suitable of cli- 
matically suitable area and has 189.43 km? of suitable wet- 
land habitat. The main parts of the wetlands are marshes, 
swamps and flood-outs. The persistence of wetlands is es- 
sential as they store atmospheric carbon in the plant roots 
and filter upstream pollutants, thus protecting coral reefs 
by holding back sediments (GIBBENS 2021, KUMAR et al. 
2017). The Agusan River is accessible by small boats and 
therefore shows high anthropogenic pressure from north- 
west to southeast portion of the marsh (GIBBENS 2021). The 
conversion of nearby areas along the river into fish ponds, 
rice fields, and/or settlements by the indigenous Manobo 
people result to areas with slightly increased anthropogenic 
pressure (Fig. 3) (Ramsar Sites Information Service 1999). 
The remaining part of the PA is exposed to low to moder- 
ate human pressure. In fact, it is even known to be the ‘least 
disturbed freshwater wetland’ in the Philippines (ASEAN 
CHM-ACB 2022). Having low human pressure and high 
suitability, the surrounding areas of Lake Mambagongon 


appear to be a hotspot for crocodiles and the lake is already 
known as a crocodile reproduction site (Tomas et al. 2009, 
VAN WEERD 2010). This PA is also a significant transit point 
for wild birds in Asia (DENR 2022) and home for 197 bird 
species as well as 53 reptile and 240 vascular plant species 
(ASEAN CHM-ACB 2022). Among the threatened species 
native to the area are the two crocodile species, the Philip- 
pine Duck (Anas luzonica), Golden-Crowned Flying Fox 
(Acerodon jubatus) and Philippine Sailfin Lizard (Hydro- 
saurus pustulatus) (DENR 2022). Philippine Crocodile 
populations were reported to occur in this PA but actually 
this revealed to be a C. porosus locality (Ross 2008). The 
coexistence of both species in the Agusan Marsh is still un- 
certain (MANALO et al. 2012). 

Legislated PAs offer a large natural habitat for the spe- 
cies with the absence or minimal presence of anthropogen- 
ic pressures. The governance of such reserves is through 
the Protected Area Management Board (PAMB) whose 
members include several stakeholders from public and 
private sectors. This management board ensures that ac- 
tivities that will be conducted within the PAs boundaries 
abide with the national laws, particularly with the NIPAS 
Act. The strict regulation of anthropogenic activities with- 
in these reserves provide a safeguard for these habitats to 
remain intact and ensure its ecological integrity. Moreover, 
legislated PAs have allocated government funds for their 
management and protection. Additionally to the top five 
reserves, we recommend Ligawasan Marsh Game Refuge 
and Bird Sanctuary on Mindanao for the establishment of 
a new reserve, which is not a declared PA yet, but in real- 
ity there is a large population of C. mindorensis present. 
The area also shows high climatic suitability and low level 
of human footprint. The declaration as a PA is highly im- 
portant for the conservation and existence of the Critically 
Endangered Philippine Crocodile. 

Most of the species occurrence records in this study 
were outside declared and legislated PAs. Establishment of 
these habitats into PAs would entail a large sum of funds 
and efforts, which usually takes years to be established. Un- 
der the Philippine Wildlife Resources Conservation and 
Protection Act of 2001, public and privately owned areas 
outside PAs which serve as a habitat for threatened species 
can be declared as a Critical Habitat. Similar to legislated 
PAs, the declaration provides a layer of protection through 
a management board which regulates activities within the 
Critical Habitat, but takes shorter time to establish. 


Crocodiles role in the ecosystem and 
human-wildlife conflicts 


Crocodiles may serve as umbrella species for their eco- 
systems, which enables the protection of other threatened 
species and entire ecosystems. In the case of the Philippine 
Crocodile, there are many other wetland-dependent spe- 
cies with threatened status that would benefit from the ex- 
pansion of existing PAs or the establishment of new ones 
in order to better protect crocodiles, including Philippine 
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natives or even endemics such as Pelochelys cantorii which 
is Critically Endangered (Brown et al. 2013), Anas luzoni- 
ca, Hydrosaurus pustulatus, Limnonectes parvus, Platyman- 
tis sierramadrensis, Sanguirana tipanan (BROWN et al. 2013, 
SANGUILA et al. 2016), which are Vulnerable, as well as di- 
verse migratory birds. 

Buco et al. (2020) and CORVERA et al. (2017) have 
shown that crocodile species native to the Philippines 
might have positive impact on the fish stocks. BRowN et al. 
(2021) suggests C. mindorensis as a potential natural pest 
control agent based on analyses of the digestive tract. In- 
vasive species such as Pomacea canaliculata (Golden Apple 
Snail) or Rattus tanezumi (Asian House Rat) are agricul- 
tural nuisances which are preyed by C. mindorensis. 

The increase of protected reserves suitable for the Phil- 
ippine Crocodile is also necessary to prevent future ex- 
tinction of the species in the wild. Unfortunately, the main 
threat for C. mindorensis concerns its habitat, particularly 
fragmentation, use and destruction. The human footprint 
index can be a useful indicator of anthropogenic expan- 
sion and habitat loss, even if it is assessed remotely and 
can slightly differ from a local scale. Our results reveal that 
many protected reserves also have larger proportions with 
moderate to high human footprint (Fig. 3). The expansion 
of agricultural land for aquaculture or for the cultivation 
of rice and sugar cane, human-settlement growth, energy 
production and lucrative mining are destroying the hab- 
itats of this species (CORVERA et al. 2017, MANALO et al. 
2018, SARMIENTO 2022, VAN WEERD & VAN DER PLOEG 
2012). Deforestation deprives them of shelter and prey re- 
sources (VAN WEERD & VAN DER PLOEG 2012), which is 
intensified by fishing activities and can have a particular- 
ly negative impact on hatchling survival (MaNALo et al. 
2015, SOMAWEERA et al. 2018). As a result, crocodiles are 
being displaced from their former habitats. They are now 
found more frequently in rice fields and near settlements, 
increasing the risk of human-crocodile conflicts (CoRVE- 
RA et al. 2017). Therefore, home ranges observed in studies 
by van WEERD et al. (2006) and VAN DE VEN et al. (2017) 
should be taken into account when selecting areas for in 
situ conservation measures and appropriate buffer zones. 
More space is necessary than is available now to deescalate 
the aggressive intraspecific, territorial behaviour in partic- 
ular of young Philippine Crocodiles (MAUGER et al. 2017, 
VAN WEERD 2010, VAN WEERD & VAN DER PLOEG 2012). 

Unlike the larger and more aggressive species Crocody- 
lus porosus, there has only been a single record of human- 
crocodile conflict in C. mindorensis in the country (CORVE- 
RA et al. 2017). For the Critically Endangered C. mindoren- 
sis, repopulating the species in the wild can be considered 
a priority conservation action. Nevertheless, regardless of 
the species, crocodiles are generally feared by most Filipi- 
nos which is a major problem for introducing crocodiles in 
suitable habitats in the country. This makes the conserva- 
tion introduction programs for the species highly complex 
involving political aspects and gathering the communities’ 
support (MANALO et al. 2015). Such complexities would 
still arise if a decision has been made to introduce the spe- 
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cies in suitable habitats within the identified PAs. On the 
other hand, the low human footprint in PAs makes them 
ideal as introduction sites with fewer human-crocodile 
interactions. In spite of such difficulties, there have been 
two conservation release programs already for the species 
(MANALO et al. 2015, VAN DE VEN et al. 2009) which proves 
that it is not impossible to introduce and repopulate the 
species in the wild. 


Acknowledgements 


We thank the Department of Environment and Natural Resourc- 
es — Biodiversity Management Bureau for providing the updated 
legislated PA maps of the Philippines. We are grateful to Morris 
FLECKS for his technical support with the figures and to VERA 
STEINBERG for commenting on a previous version of the manu- 
script. Furthermore, we would like to thank ARNE SCHULZE and 
Mona VAN SCHINGEN-KuaN for editorial support, reviewing the 
manuscript and providing advice. 


References 


ADAMS, W. M., R. AVELING, D. BROCKINGTON, B. Dickson, J. EL- 
LIOTT, J. HUTTON, D. Rog, B. VIRA & W. WOLMER (2004): Bio- 
diversity conservation and the eradication of poverty. - Sci- 
ence, 306: 1146-1149. 


AKMAD, M. & C. POMARES (2008): Wild Crocodylus mindorensis 
at Tambad, Carmen, North Cotabato. - Unpublished BSc (Ag- 
riculture) thesis. Kabacan, North Cotabato, Philippines, De- 
partment of Animal Science, College of Agriculture, Univer- 
sity of Southern Mindanao. 


ALLOUCHE, O., A. TSoAR & R. KADMON (2006): Assessing the ac- 
curacy of species distribution models: prevalence, kappa and 
the true skill statistic (TSS). - Journal of Applied Ecology, 43: 
1223-1232. 

AMBAL, R., M. Duya, M. Cruz, O. Coroza, S. VERGARA, N. DE 
SILVA, N. MOLINYAWE & B. TABARANZA (2012): Key Biodiver- 
sity Areas in the Philippines: Priorities for Conservation. - 
Journal of Threatened Taxa, 4: 2788-2796. 


ASEAN CHM-ACB (2022): Agusan Marsh Wildlife Sanctuary. - 
https://asean.chm-cbd.net/agusan-marsh-wildlife-sanctuary, 
accessed 16 February 2022. 


Brnapay, J. W. B., R. I. MANALO, P. C. BALTAZAR & F. P. MAGAL- 
LANES (2020): Crocodylus mindorensis: Habitat use. - Herpe- 
tological Review, 51: 319-321. 


Brnapay, J. W. B., R. I. MANALO, M. A. Braaalis, L. K. M. PALAO, 
L. R. P. Mojica, O. B. S. NICOPIOR & P. C. BALTAZAR (2020): 
Determining habitat suitability for the identification of prior- 
ity conservation sites for Indo-Pacific crocodile (Crocodylus 
porosus Schneider, 1801) in the Philippines. - Philippine Jour- 
nal of Science, 150: 333-344. 


Brown, J. C., M. H. SHIRLEY, A. YoG-yoGc, M. VAN WEERD, M. 
G. BALBAs, B. A. TARUN & C. D. SILER (2021): Use of diet and 
body condition assessments as intermediate indicators of 
translocation success in the Critically Endangered Philippine 
crocodile (Crocodylus mindorensis). - Aquatic Conservation: 
Marine and Freshwater Ecosystems, 31: 2817-2829. 


Brown, R. M., C. D. SILER, C. H. OLIVEROS, L. J. WELTON, A. 
Rock, J. SWAB, M. vAN WEERD, J. V. BEIJNEN, E. Jose, D. RoD- 


Suitable protected areas for Crocodylus mindorensis 


RIGUEZ, E. Jose & A. C. DiEsmos (2013): The amphibians and 
reptiles of Luzon Island, Philippines, VIII: the herpetofauna of 
Cagayan and Isabela Provinces, northern Sierra Madre moun- 
tain range. - ZooKeys, 266: 1-120. 


Bucor, A. A., R. I. MANALO, A. C. ALCALA & P. S. ASPILLA 
(2020): Crocodilians and fisheries in the Philippines: revis- 
iting Fittkau’s hypothesis. - Philippine Journal of Systematic 
Biology, 14: 9. 


CBD (2018): COP14 - Breakfast: Ministerial Meeting of Like- 
minded Megadiverse Countries (LMMC) Group (15 Novem- 
ber). - https://www.cbd.int/conferences/2018/parallel-meet- 
ings/like-minded-megadiverse-countries-Immc, accessed 13 
September 2022. 


CEPF (2001): Philippines. - https://www.cepf.net/our-work/bi- 
odiversity-hotspots/philippines, accessed 13 September 2022. 


Corp, A. & D. RODDER (2011): Inclusion of habitat availability in 
species distribution models through multi-temporal remote- 
sensing data? - Ecological Applications, 21: 3285-3298. 


CorverA, M., R. MANALO & M. Aquino (2017): People and croc- 
odiles sharing one environment: an analysis of local human- 
crocodile conflict management strategies in the Philippines. 
— Journal of Animal Science and Research, 1: 1-6. 


De Marco, P. & C. C. NÓBREGA (2018): Evaluating collinearity 
effects on species distribution models: An approach based on 
virtual species simulation. - PloS One, 13: e0202403. 

DENR-BMB (2020): Statistics. - https://bmb.gov.ph/index.php/ 

resources/statistics, accessed 5 September 2021. 

DENR-BMB (2021): Republic Acts. - Department of Environ- 

ment Natural Resources, Biodiversity Management Bureau, 

Philippines. 

DENR (2015): B+ WISER - Program Site Northern Sierra Madre 

Natural Park. - https://forestry.denr.gov.ph/b+wiser/index. 

php/sites/nsmnp, accessed 15 February 2021. 

DENR (2022): Agusan Marsh Wildlife Sanctuary. - http://den- 
rpenroads.com/index.php/priority-programs/national-green- 
ing-program/protected-area-biodiversity-and-ecotourism, 
accessed 1 March 2022. 

DENR (2023): Lake Lanao Watershed. — https://forestry.denr.gov. 
ph/inremp/lanao.html, accessed 22 February 2023. 

ELITH, J. & C. H. GRAHAM (2009): Do they? How do they? WHY 
do they differ? On finding reasons for differing performances 
of species distribution models. - Ecography, 32: 66-77. 


ELITH, J., C. H. GRAHAM, R. P. ANDERSON, M. Dupik, S. FERRIER, 
A. GUISAN, R. J. HIJMANS, F HUETTMANN, J. R. LEATHWICK 
& A. LEHMANN (2006): Novel methods improve prediction of 
species’ distributions from occurrence data. - Ecography, 29: 
129-151. 


ELITH, J., M. KEARNEY & S. PHILLIPS (2010): The art of modelling 
range-shifting species. - Methods in ecology and evolution, 
1: 330-342. 

Fors, M., A. CUENA-LOMBRANA, G. FENU & G. BACCHETTA 
(2018): Using species distribution models at local scale to 
guide the search of poorly known species: Review, methodo- 
logical issues and future directions. - Ecological Modelling, 
385: 124-132. 


FPE (2013): Biodiversity - Hub of Life: Species Diversity in the 
Philippines. - https://fpe.ph/biodiversity.html/view/hub-of- 
life-species-diversity-in-the-philippines, accessed 13 Septem- 
ber 2021. 


GIBBENS, S. (2021): The world’s wetlands are slipping away. This 
vibrant sanctuary underscores the stakes. - National Geo- 
graphic - Environment. - https://www.nationalgeograph- 
ic.com/environment/article/world-wetlands-are-slipping- 
away-agusan-marsh-underscores-stakes, accessed 5 October 
2021. 


GinaL, P, W. C. Tan & D. RÖDDER (2022): Invasive risk assess- 
ment and expansion of the realized niche of the Oriental Gar- 
den Lizard Calotes versicolor species complex (Daudin, 1802). 
- Frontiers of Biogeography, 14: e54299 


GRDC (2020): WMO Basins and Sub-Basins / Global Runoff - 
Data Centre GRDC. Koblenz, Germany. 


GUMBRICHT, T., R. M. RoMAN-CuEsTA, L. VERCHOT, M. HEROLD, 
F. WITTMANN, E. HOUSEHOLDER, N. HEROLD & D. MURDI- 
YARSO (2017): An expert system model for mapping tropical 
wetlands and peatlands reveals South America as the largest 
contributor. - Global Change Biology, 23: 3581-3599. 


HIJMANS, R., S. CAMERON, J. PARRA, P. JONES & A. JARVIS (2005): 
WorldClim - Very High Resolution Interpolated Climate Sur- 
faces for Global Land Areas. Version 1.4 (release 3). — Muse- 
um of Vertebrate Zoology - University of California, Berkeley, 
www. worldclim.org. 


InLow, F, R. Bonke, T. HARTMANN, P. GEISSLER, N. BEHLER & 
D. RODDER (2015): Habitat suitability, coverage by protected 
areas and population connectivity for the Siamese crocodile 
Crocodylus siamensis Schneider, 1801. - Aquatic Conservation: 
Marine and Freshwater Ecosystems, 25: 544-554. 


IUCN (2021a): The IUCN Red List of Threatened Species. Ver- 
sion 2021-2. — https://www.iucnredlist.org, accessed 5 Septem- 
ber 2021. 


IUCN (2021b): Protected Area Categories. - https://portals.iucn. 
org/library/sites/library/files/documents/pag-021.pdf, ac- 
cessed 5 August 2021. 


JAISANKAR, I., A. VELMURUGAN & C. SIVAPERUMAN (2018): Bio- 
diversity conservation: issues and strategies for the tropical is- 
lands. - pp. 525-552 in: SIVAPERUMAN, C., A. VELMURUGAN, 
A. K. SINGH & I . JAISANKAR (eds): Biodiversity and climate 
change adaptation in tropical islands. - Elsevier. 


KUMAR, R, R. J. McCINNEs, M. EVERARD, R. C. GARNER, K. A. A. 
KULINDWA, H. WITTMER & D. INFANTE Mata (2017): Integrat- 
ing multiple wetland values into decision-making. - Ramsar 
Policy Brief No. 2., Gland, Switzerland. 


Loso, J. M., A. JIMÉNEZ-VALVERDE & R. REAL (2008): AUC: a 
misleading measure of the performance of predictive distribu- 
tion models. - Global Ecology and Biogeography, 17: 145-151. 


MANALO, R. (2007): Occurrence of Crocodylus mindorensis in the 
Cordillera Central. Abra Province, Luzon Island. - National 
Museum Papers, 14: 109-115. 


MANALO, R., A. C. ALCALA & A. BucoL, W. T. BELO, V. P. MER- 
CADO & P. C. BALTAZAR (2015): Philippine Crocodile Release 
Program: An update on the status of the introduced Philippine 
Crocodiles (Crocodylus mindorensis) in Paghungawan Marsh, 
Siargao Island Protected Landscapes and Seascapes, Pilar, 
Surigao Del Norte, Philippines. - First East & South East Re- 
gional Meeting of the IUCN-SSC Crocodile Specialist Group, 
Siem Reap, Cambodia, 25-29 May, 2015. 


MANALO, R, J. W. Binapay, P. BALTAZAR, M. CORVERA & C. J. 
LADIANA (2019): New population record of Philippine Croco- 
dile (Crocodylus mindorensis Schmidt, 1935) in southern Phil- 
ippines. - Crocodile Specialist Group Newsletter, 38: 7-10. 


39 


SABINE HARRER et al. 


MANALO, R. I., W. T. BELO, V. P. MERCADO, B. O. SOLCO & J. A. J. 
BINAN (2012): Distribution and status of crocodiles in Agusan 
Marsh, eastern Mindanao, Philippines. - pp. 50-57 in: Pro- 
ceeding of the 21* Working Meeting of the IUCN-SSC Croco- 
dile Specialist Group (org.): Crocodiles. - IUCN, Gland, Swit- 
zerland. 


MANnaLo, R. I., C. C. POMARES, V. P. MERCADO, W. T. BELO, G. 
SALJAY & T. D. Turas (2013): A new distribution record for the 
Philippine crocodile (Crocodylus mindorensis, Schmidt 1935). 
- World crocodile conference. Proceedings of the 22"! work- 
ing meeting of the IUCN-SSC Crocodile Specialist Group, Ne- 
gombo, Sri Lanka, 21-23 May 2013. 


Mana_o, R. I., E. A. TABAYAG & P. C. BALTAZAR (2018): Con- 
servation milestones of the Critically Endangered Philippine 
Crocodile (Crocodylus mindorensis, Schmidt 1935). - Sylvat- 
rop, The Technical Journal of Philippine Ecosystems and Nat- 
ural Resources, 28: 31-38. 


MAUGER, L. A., E. VELEZ, M. S. CHERKISS, M. L. BRIEN, F. J. Maz- 
ZOTTI & J. R. SPOTILA (2017): Conservation genetics of Ameri- 
can crocodile, Crocodylus acutus, populations in Pacific Costa 
Rica. - Nature Conservation, 17: 1. 


MAazzorTTI, F, G. BEST, L. BRANDT, M. CHERKISS, B. JEFFERY & 
K. Rice (2009): Alligators and crocodiles as indicators for res- 
toration of Everglades ecosystems. - Ecological Indicators, 9 
(Supplement): 137-149. 


MCALLISTER, D. E., J. E Craic, N. Davipson, S. DELANY & 
M. SEDDON (2001): Biodiversity Impacts of Large Dams. - 
Background paper Nr 1, Prepared for IUCN/UNEP/WCD. - 
IUCN/UNEP. 


Mcteop, D. S., C. D. SILER, A. C. Diesmos, M. L. DiEsmos, V. S. 
GARCIA, A. O. ARKONCEO, K. L. BALAQuitT, C. C. Uy, M. M. 
VILLASERAN, E. C. YARRA & R. M. Brown (2011): Amphibi- 
ans and reptiles of Luzon Island, V: The herpetofauna of Angat 
Dam Watershed, Bulacan Province, Luzon Island, Philippines. 
- Asian Herpetological Research, 2: 177-198. 


Merrow, C., M. J. SMITH & J. A. SILANDER JR (2013): A practi- 
cal guide to MaxEnt for modeling species’ distributions: what 
it does, and why inputs and settings matter. - Ecography, 36: 
1058-1069. 


Mu, Q., E A. HEINSCH, M. ZHAO & S. W. RUNNING (2007): De- 
velopment of a global evapotranspiration algorithm based on 
MODIS and global meteorology data. - Remote sensing of En- 
vironment, 111: 519-536. 


Namı, B., N. A. HAMM, T. A. GROEN, A. K. SKIDMORE & A. G. 
ToxopEus (2014): Where is positional uncertainty a problem 
for species distribution modelling? - Ecography, 37: 191-203. 


PEREZ, P. L. & BUKLURAN - The Philippine ICCA Consortium. 
(2018): Living with the problem of national parks: Indigenous 
critique of Philippine environmental policy. - Thesis Eleven, 
145: 58-76. 

PERIA, E. V. (2014): The Philippines - Monograph on Stakeholder 
Participation in the NBSAP Revision Process. 


PHILLIPS, S., M. Dupik & R. E. SCHAPIRE (2016): Maxent soft- 
ware for modeling species niches and distributions (Version 
3-4.0). - https://biodiversityinformatics.amnh.org/open_ 
source/maxent/. 


PHILLIPS, S. J., R. P. ANDERSON & R. E. SCHAPIRE (2006): Maxi- 
mum entropy modeling of species geographic distributions. - 
Ecological Modelling, 190: 231-259. 


40 


PHILLIPS, S. J. & M. Dupfxk (2008): Modeling of species distribu- 
tions with Maxent: new extensions and a comprehensive eval- 
uation. - Ecography, 31: 161-175. 


PHILLIPS, S. J., M. DUDÍK, J. ELITH, C. H. GRAHAM, A. LEHMANN, 
J. LEATHWICK & S. FERRIER (2009): Sample selection bias and 
presence-only distribution models: implications for back- 
ground and pseudo-absence data. - Ecological Applications, 
19: 181-197. 


PSA (2019): Philippine Statistics Authority. - https://psa.gov.ph/ 
article/terms-use, accessed 5 October 2021. 


QGIS.org (2021): QGIS Geographic Information System, QGIS 
Association. — http://www.qgis.org. 


R Core Team (2020): R: A language and environment for statisti- 
cal computing. - R Foundation for Statistical Computing, Vi- 
enna, Austria. 


Ramsar Sites Information Service (1999): Agusan Marsh Wildlife 
Sanctuary. — https://rsis.ramsar.org/ris/1009, accessed 5 Octo- 
ber 2021. 


RODDER, D., J. O. ENGLER, R. BONKE, F. WEINSHEIMER & W. PER- 
TEL (2010): Fading of the last giants: an assessment of habitat 
availability of the Sunda gharial Tomistoma schlegelii and cov- 
erage with protected areas. - Aquatic Conservation: Marine 
and Freshwater Ecosystems, 20: 678-684. 


Ross, C. (1982): Philippine crocodile project. Final report. - 
Smithsonian Institute and World Wildlife Fund, Washington. 


Ross, C. A. (ed.) (2008): Proceedings of the forum on crocodiles 
in the Philippines. - National Museum Papers, 14. National 
Museum of the Philippines, Manila. 


Ross, C. A. & A. C. ALCALA (1983): Distribution and status of the 
Philippine crocodile (Crocodylus mindorensis). - Kalikasan, 
12: 169-173. 


SANGUILA, M. B., K. A. Coss, C. D. SILER, A. C. Diesmos, A. C. 
ALCALA & R. M. Brown (2016): The amphibians and reptiles 
of Mindanao Island, southern Philippines, II: the herpetofau- 
na of northeast Mindanao and adjacent islands. - ZooKeys, 
624: 1-132. 


SARMIENTO, B. S. (2022): Philippine groups slam ‘cruel Christmas 
gift’ as open-pit mining ban is lifted. - Mongabay Series: Land 
rights and extractives, https://news.mongabay.com/2022/01/ 
philippine-groups-slam-cruel-christmas-gift-as-open-pit-mi- 
ning-ban-is-lifted/, accessed 11.06.2023. 


SCHARLEMANN, J. P. W., D. Benz, S. I. Hay, B. V. Purses, A. J. Ta- 
TEM, G. R. W. WINT & D. J. RoGERS (2008): Global data for 
ecology and epidemiology: a novel algorithm for temporal 
Fourier processing MODIS data. - PloS One, 3: e1408. 


SHABANI, F., L. KUMAR & M. AHMADI (2018): Assessing accuracy 
methods of species distribution models: AUC, specificity, sen- 
sitivity and the true skill statistic. - Global Journal of Human 
Social Science, 18: 6-18. 


SOMAWEERA, R., M. L. BRIEN, S. G. PLATT, C. MANOLIS & B. L. 
WEBBER (2018): Direct and indirect interactions with vegeta- 
tion shape crocodylian ecology at multiple scales. - Freshwa- 
ter Biology, 64: 257-268. 


Tan, W. C., P. GINAL, A. G. RHODIN, J. B. IVERSON & D. RODDER 
(2022): A present and future assessment of the effectiveness 
of existing reserves in preserving three critically endangered 
freshwater turtles in Southeast Asia and South Asia. - Fron- 
tiers of Biogeography, 14: e50928. 


Suitable protected areas for Crocodylus mindorensis 


Tomas, R. C., J. B. MANUTA & V. D. Rosa (2009): Women, Wa- 
ter and the Marsh: Adaptation Pathways of Agusan Marsh In- 
digenous Communities in Mindanao, Southern Philippines. - 
JPAIR Multidisciplinary Research Journal, 2: 1. 


TsuyaMa, I., M. Hica, K. Nakao, T. Matsui, M. Horikawa & 
N. TANAKA (2015): How will subalpine conifer distributions 
be affected by climate change? Impact assessment for spatial 
conservation planning. - Regional Environmental Change, 15: 
393-404. 

UETZ, P., P. FREED, R. AGUILAR, F. REYES & J. Hošek (eds) (2023): 
The Reptile Database. - http://www.reptile-database.org, ac- 
cessed 05.01.2023. 


UNDP (2021): About the Philippines. - https://www.ph.undp. 
org/content/philippines/en/home/countryinfo/, accessed 5 
September 2021. 


UNEP-WCMC & IUCN (2021): Protected Planet: The World Da- 
tabase on Protected Areas (WDPA), March 2021. - UNEP- 
WCMC and IUCN, Cambridge, UK. 


USAID (2021): The Philippines - Environment. - https://www. 
usaid.gov/philippines/energy-and-environment, accessed 5 
September 2021. 


VAN DE VEN, A., J. P. GUERRERO, D. G. RODRIGUEZ, S. P. TELAN, 
M. G. Barsas, B. A. TARUN, M. v. WEERD, J. PLOEG, Z. WIJT- 
TEN & F LINDEYER (2009): Effectiveness of head starting to 
bolster Philippine crocodile Crocodylus mindorensis. - Con- 
servation Evidence, 6: 111-116. 


VAN DE VEN, W. S. TELAN, E. José, F. LINDEYER, N. TUBBS, J. 
DE JonGe, D. RODRIGUEZ, M. BALBAS, J. GUERRERO, M. VAN 
WEERD, J. VAN DER PLOEG & H. DE IONGH (2017): Movements 
and home ranges of Philippine crocodiles in San Mariano, 
lsabela. - NMP Journal of Natural History, 2: 38-45. 


VAN DER PLOEG, J., M. VAN WEERD, A. B. MASIPIQUENA & G. A. 
Persoon (2011): Illegal logging in the Northern Sierra Madre 
Natural Park, the Philippines. - Conservation and Society, 9: 
202-215. 

vAN WEERD, M. (2010): Philippine Crocodile Crocodylus mindo- 
rensis. - pp. 71-78 in: MANOLIs, S. C. & C. STEVENSON (eds): 
Crocodiles: Status survey and Conservation Action Plan. 
Third edition. - Crocodile Specialist Group, Darwin. 


VAN WEERD, M., C. POMARO, J. DE LEON, R. ANTOLIN & V. MER- 
CADO (2016): Crocodylus mindorensis. - The IUCN Red List 
of Threatened Species 2016: e.T5672A3048281, https://dx.doi. 
org/10.2305/IUCN.UK.2016-3.RLTS.T5672A3048281.en, ac- 
cessed 13 September 2021. 


VAN WEERD, M. & J. VAN DER PLOEG (2012): The Philippine croco- 
dile: ecology, culture and conservation. - Mabuwaya Founda- 
tion, Cabagan, Isabela, Philippines. 


VAN WEERD, M., J. VAN DER PLOEG, D. RODRIGUEZ, J. GUERRERO, 
B. Tarun, S. TELAN & J. DE JONGE (2006): Philippine croco- 
dile conservation in Northeast Luzon: an update of population 
status and new insights into Crocodylus mindorensis ecology. 
- The 18" Working Meeting of the IUCN-SSC Crocodile Spe- 
cialist Group, Montelimar, France. IUCN, Gland, Switzerland 
and Cambridge, UK. 


VENTER, O., E. W. SANDERSON, A. MAGRACH, J. R. ALLAN, J. 
BEHER, K. R. JONES, H. P. PossINGHAM, W. F. LAURANCE, P. 
Woon, B. M. FEKETE, M. A. LEvy & J. E. WATSON (2018): Last 
of the Wild Project, Version 3 (LWP-3): 2009 Human Foot- 
print, 2018 Release. - NASA Socioeconomic Data and Appli- 
cations Center (SEDAC), Palisades, NY, USA. 


VON RINTELEN, K., E. ARIDA & C. HAUSER (2017): A review of 
biodiversity-related issues and challenges in megadiverse In- 
donesia and other Southeast Asian countries. — Research Ideas 
and Outcomes, 3: e20860. 


WARREN, D. L. & S. N. SEIFERT (2011): Ecological niche mode- 
ling in Maxent: the importance of model complexity and the 
performance of model selection criteria. - Ecological Applica- 
tions, 21: 335-342. 


Supplementary data 


The following data are available online: 


Supplementary Table S1. Final ranking containing detailed infor- 
mation on all 248 PAs including wetland categories, lower rank- 
ings, and further statistical values. 


Supplementary Table S2. List of pre-processed variables obtained 
from WorldClim and MODIS before calculating multi-colline- 
arity. 


41 


